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Abstract 


A  phenomenological  theory  is  developed  for  the 
propagation  of  plane  electromagnetic  waves  in  a  deformed 
non-absorbing  centrosymmetric  isotropic  material.  It  is 
assumed  that  the  dielectric  constant  and  specific  reluctance 
matrices  depend  on  the  deformation  gradients  at  the  instant 
of  measurement.  The  theory  is  formulated  from  both  the 
Eulerian  and  Lagrangian  standpoints. 
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1 .  Introduction 


In  this  paper  we  consider  the  propagation  of  plane 
electromagnetic  waves  in  a  non-absorbing  material  which  is 
subjected  to  finite  deformations.  It  is  assumed  that  the 
material  is  isotropic  when  undeformed  and  when  no  electro¬ 
magnetic  fields  are  present  and  that  it  is  centrosymmetric. 

The  theory  is  formulated  from  both  Eulerian  and  Lagrangian 
points-of-view.  The  latter  formulation  rests  on  the  Lagrangian 
formulation  of  Maxwell's  equations  for  a  deformed  material  due 
to  Walker,  Pipkin  and  Rivlin  [1] . 

In  each  case  the  assumption  is  made  that  the  material 
is  linear  with  respect  to  electromagnetic  effects,  but  that 
the  dielectric  constant  and  specific  reluctance  matrices  may 
depend  on  the  displacement  gradients  in  the  material.  It 
follows  from  the  isotropic  character  of  the  material  that  the 
dielectric  constant  and  specific  reluctance  matrices  are 
isotropic  matrix  functions  of  the  Cauchy  and  Finger  strains, 
accordingly  as  the  Lagrangian  or  Eulerian  formulation  is  adopted 
and  may  be  expressed  in  terms  of  these  in  canonical  forms.  In 
each  case  we  obtain  from  the  constitutive  equations  and  Maxwell's 
equations  a  secular  equation  for  the  determination  of  the  slow¬ 
ness  of  a  plane  electromagnetic  wave,  propagating  in  an  arbitrary 
direction  in  a  material  which  is  subjected  to  a  pure  homogeneous 
deformation. 

We  pursue  the  study  of  this  equation  in  the  Eulerian 
case  and  obtain  the  six  principal  slownesses.  It  is  found 
that  there  is  a  relation  between  these  six  slownesses.  In  the 
case  when  only  the  dielectric  constant  or  only  the  specific 
reluctance  depends  on  the  deformation,  this  single  relation  is 
replaced  by  three  relations. 
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We  then  discuss  the  propagation  of  the  electromagnetic 
wave  in  any  direction  in  a  principal  plane.  In  §§  5  and  6  we 
consider  propagation  in  a  material  which  is  subjected  to  shear¬ 
ing  deformations. 

Finally  in  §  7  we  consider  the  application  of  the 
theory  to  materials  for  which  the  dielectric  constant  and 
specific  reluctance  matrices  depend  on  the  history  of  the  de¬ 
formation,  but  in  which  the  deformation  is  held  constant. 
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2.  The  constitutive  equations 
(a)  Eulerian  formulation 


We  consider  a  body  to  undergo  a  deformation  which 
is  described  in  a  rectangular  cartesian  coordinate  system  x 
by 


xi  =  xi(t)  =  xi(XA/t)/  (2.1) 


where  x.^  is  the  position  in  the  system  x,  at  time  t,  of  a 

particle  which  was  at  in  the  same  system  at  a  reference 

time  t  . 
o 

We  make  the  constitutive  assumption  that  the  electric 
displacement  field  d^  at  time  t,  depends  only  on  the  electric 
field  e^  and  deformation  gradients  x^  A,  measured  at  the  particle 
considered  at  time  t.  We  also  assume  that  the  dependence  of 
d.  on  e  is  linear.  We  make  the  analogous  constitutive  assumption 

i  p  “ 

that  the  magnetic  induction  field  b^,  at  time  t,  depends  only  on 
the  magnetic  intensity  field  hp  and  deformation  gradients  A, 
the  dependence  on  the  former  being  linear. 

If  the  material  is  isotropic  in  its  reference  state, 
it  follows  [2]  that 


d.  =  k.,e,  and  h.  =  w.,b.f  (2.2) 

i  ij  j  i  ij  y 


where  k  ,  the  dielectric  constant  tensor,  and  u>.  ,  the  specific 

1  J  1  J 

reluctance  tensor,  are  given  by 


k,  ,  =  k_6_.  J  +  k,c_,  ,  +  k_c.,.c,. 


i  j 


°  ij 


1  ij 


2  ik  kj 


and 


(2.3) 


(D.  ,  =  W  ($.,  +  (O.c,  ,  +  W0C..C,  ., 
ij  o  ij  1  ij  2  ik  kj' 
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where  c.^  is  the  Finger  strain  tensor  defined  by 


c.  .  =  x.  .x,  .  -  6 .  . . 
lj  1 ,  A  j  ,A  lj 


(2.4) 


In  (2.3),  k  ,  k.,  kn,  03,0).,  to-  are  functions  of  the  invar- 
o  1  2  o  1  2 


2  3 

iants  tr  c,  tr  c  ,  tr  c  ,  where  c  = 


the  notation  k  = 


k. 

1 1 


Cij 


Introducing 
,  e  =  (e^)  ,  with  analogous  meanings 


for  other  bold-face  symbols,  we  may  rewrite  (2.2)  as 


d  =  k.e  and  h  =  to.b, 

where 

k  =  k  I  +  knc  +  k0c2, 

~  o~  1~  2~ 

2 

(0  =  (i)  I  +  ID.C  +  to0C  , 

~  o~  1 ~  2 ~ 

and  I  denotes  the  unit  matrix. 

For  a  plane  electromagnetic  wave,  adopting  the  usual 
complex  notation,  we  may  write  e,  h,  d,  b  in  the  form 

(e,  h,  d,  b)  =  (e,  h,  d,  b)  e  (2.7) 

<y  *w  ^  A,  A#  ^  ^  A, 

where  e,  h,  d,  b  and  s  are  vectors  which  may  be  real,  imaginary, 
or  complex  constants.^  s  is  the  complex  slowness  of  the  wave 
and  a)  is  its  angular  frequency.  Then,  the  constitutive  equations 
(2.5)  become 

d  =  k.e,  h  =  to.b.  (2.8) 


(2.5) 


(2.6) 


1*  . 

We  will  see  later  that  for  the  constitutive  equation  discussed 
here,  the  case  when  s  is  complex  can  be  ruled  out. 


(b)  Lagranqian  formulation 


An  alternative  formulation  may  be  attained  in  the 
following  way.  In  accord  with  Walker,  Pipkin  and  Rivlin 
[1] ,  we  define  the  reduced  fields  E,  H,  B,  D  by  the  equations 


E  =  F*e,  H  =  F?h, 

B  =  (det  F)  F-1.b  and  D  =  (det  F)  F-1.d, 


(2.9) 


where  the  notation 


!l 


(2.10) 


is  used  and  the  star  denotes  the  transpose.  The  constitutive 
assumptions  made  as  a  basis  for  the  Eulerian  formulation  are 
equivalent  to  the  assumptions  that  D  and  B  are  linear  functions 

*v 

of  E  and  H  respectively  and  both  D  and  B  depend  on  F. 

Then,  the  assumption  that  the  material  is  isotropic 
in  its  reference  state  leads  to  the  conclusion  that 


D  =  K.E  and  H  =  ft.B 

and  K  and  ft  are  expressible  in  the  forms 

K=K  I  +K.C+  K0C2 
~  o~  1~  2~ 


(2.11) 


and 


ft  =  ft  I  +  ft.  C 
o~  1~ 


ft2c2, 


(2.12) 


where  C  is  the  Cauchy  strain  defined  by 


-  JJ  ABjl  |  i,A  i  ,B  AB 

2 

and  Kq,  ,  Kg,  fiQ,  fi.^,  fig  are  functions  of  tr  C,  tr  C  ,  tr 
The  relations  between  K  ,  fi  (a=0,l,2)  and  k  ,  to 

0*  Ot  OC  OC 

(01=0,1,2)  can  be  derived ,  However,  the  algebra  involved  is 
somewhat  cumbersome. 

Now,  we  consider  the  electromagnetic  wave  for  which, 
adopting  the  usual  complex  notation, 

(E,  H,  D,  B)  =  (E,  H,  D,  3)elW^,£“t^,  (2.14) 

where  E,  H,  D,  B  and  S  may  be  real,  imaginary  or  complex 
constant  vectors.  We  obtain  from  (2.11) 

D  =  K.E ,  H  =  fi.B.  (2.15) 

We  note  that  if  the  electromagnetic  fields  e,  h,  d,  b  corres¬ 
pond  to  a  plane  wave,  i.e.,  are  of  the  form  (2.7),  the  derived 
electromagnetic  fields  E,  H,  D,  B  will  not,  in  general,  have 
the  form  (2.14). 


t  o 


-  /  - 


3.  Derivation  of  the  secular  equation 


(a)  Eulerian  formulation 


Maxwell's  equations  may  be  written  in  the  form 


curl  e  =  -3b/3t,  curl  h  =  3d/3t, 


(3.1) 


where 


(curl  e).  =  e...e.  , 

~  i  i jk  k,j 


(3.2) 


Introducing  (2.7),  we  obtain 


eijksjek  V  £ijksjhk  “di* 


(3.3) 


Eliminating  e,  h  and  b  from  equations  (3.3)  and  (2.8),  we 
obtain , 


6 .  .  +  e .  e  s  s  w  (k-1)  . 

lj  ipq  mrs  p  r  qm  ~  sj 


dj  =  °. 


(3.4) 


Alternatively  eliminating  e,  h  and  d  from  equations  (3.3) 
and  (2.8)  we  obtain 


<5 .  .  +  s .  s  s  s  (k  1)  w  , 
lj  ipq  mrs  p  r  ~  qm  sj 


b  =  0. 

u 


(3.5) 


Equation  (3.4)  yields  a  non- trivial  solution  for 
d  and  (3.5)  yields  a  non-trivial  solution  for  b  provided 
that 


k. ,  +  e .  e.  ssw 

ipq  jsr  p  r  qs 


0. 


(3.6) 
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are 


For  the  wave  (2.7),  the  planes  of  constant  phase 


s  .x  =  constant 


(3.7) 


and  the  planes  of  constant  amplitude  are 


s  .x  =  constant. 


(3.8) 


In  the  particular  case  when  these  are  the  same,  we  may  write 


s  =  sn. 


(3.9) 


where  n  is  a  (real)  unit  vector  and  s  is  a  constant  which  may 
be  real,  imaginary,  or  complex.  Then,  equations  (3.4)  and  (3.5) 
become 


and 


O  _  3 

<5 .  .  +  s“e .  e  n  n  a)  (k  )  . 
lj  ipq  rars  p  r  qm  ~  sj 


2  -1 
6..  +  se.  e  nn(k  )  w. 
ij  ipq  mrs  p  r  -  qra  sj 


dj =  0 


(3.10) 


=0, 


and  (3.6)  yields  the  secular  equation  for  the  complex 
slowness  s. 


k.  .  +  s  e.  e .  n  n  w  =  0 . 
1  i j  ipq  jsr  p  r  qs 1 


(3.11) 


We  shall  call  the  direction  of  n  the  direction  of  propagation 
of  the  wave. 

It  is  shown  in  the  Appendix  that  (3.11)  may  be  written 
as 


k  2 

4>s  -  tps  +9  =  0, 


(3.12) 


where 


<J>  =  (n.k.n)  (n.oT^.n)  det  co, 

^  =  n.{(tr  k  w)k  -  k  to  k}.n,  (3.13) 

6  =  det  k. 


Equations  (3.10)  may  be  simplified  slightly  by 
choosing  the  reference  system  so  that  the  unit  normal  to  the 
wave-front  is  in  the  direction  of  the  x^-axis,  i.e.,  so  that 
ni  =  ^i3*  Equations  (3.3)  then  yield,  with  (3.9), 

d3  =  b3  =  0  (3.14) 


and  equations  (3.10)  become 


and 


(3.15) 


s^eayepT^~  ^  yTWp$ 


bB  0  ' 


where  Greek  indices  take  the  values  1,  2  and  e  0  denotes  the 

ccp 

two-dimensional  alternating  symbol. 

It  is  evident  from  (3.14)  and  (3.15)  that  the  wave 
is,  in  general,  polarized  elliptically  with  its  electric 
displacement  and  magnetic  induction  fields  in  planes  normal 
to  the  direction  of  propagation.  It  then  follows  from  (2.8) 
that  e  and  h  are,  in  general,  not  perpendicular  to  the 
direction  of  the  propagation. 
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p 

If  both  of  the  values  of  s  given  by  (3.17)  are  positive, 

then  we  obtain  two  positive  values  of  s  and  two  negative 

values.  This  corresponds  to  the  possibility  of  two  waves 

in  the  positive  direction  of  n  and  two  waves  in  the  negative 

direction.  If,  on  the  other  hand,  <p,ip  and  0  are  such  that 

2 

for  any  n,  one  of  the  values  of  s  ,  given  by  (3.17)  ,  is 
negative,  the  corresponding  values  of  s  are  imaginary.  The 
material  would  then  be  inherently  electromagnetically  unstable 
in  the  state  of  deformation  considered. 
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(b)  Lagrangian  formulation 

It  has  been  pointed  out  by  Walker,  Pipkin  and 
Pdvlin  [1]  that  in  terms  of  the  derived  electromagnetic 
fields  E,  H,  D  and  B,  Maxwell's  equations  can  be  written  as 

Curl  E  =  -  9B/9t,  Curl  H  =  9D/3t,  (3.19) 


where 


(Curl  E)a  eA3CEc1B‘ 


(3.20) 


Introducing  (2.14)  into  (3.19) ,  we  obtain 


£abcsbec  ba'  £abcsbhc  ~da* 


(3.21) 


Eliminating  E,  H  and  B  from  (3.21)  and  (2.15),  we  obtain 


6AB  +  eAPQeMRSSPSRfiQM(?  ^  SB 


db  =  0. 


(3.22) 


Again,  eliminating  E,  H  and  D  from  (3.21)  and  (2.15)  we 
obtain 


6AB  +  eAPQeMRSSPSR }QMQSB 


BB  =  °* 


(3.23) 


Again,  if  the  planes  of  constant  amplitude  and  phase 
in  the  X-space  are  the  same  and  N  is  the  unit  normal  perpendicular 
to  this  plane,  we  may  write  analogously  with  (3.9), 


S  -■  SN, 


(3.24) 


where  S  is  a  constant  which  may  be  real,  imaginary,  or 
complex.  Then,  equations  (3.22)  and  (3.23)  become 
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6AB  +  S  £APQ£MRSNPNRfiQM(?  *  SB 


Dfi  =  0 


and 


(3.25) 


6AB  +  S  £APQ£MRSNPNR(~  *  QM^SB 


The  secular  equation  for  S  is 


BB  =  0, 


lKAB  +  s  £APQ£BMLNPNLQQmI  °* 


(3.26) 


Following  a  procedure  similar  to  that  used  in  the  Appendix 
to  derive  (3.12),  we  can  express  (3.22)  in  the  form 


4>S 


h 


TS2  +  0  =  0, 


(3.27) 


where 


$  = 


(N.K,N)  (N.Q~1.N)det  Q, 

M  <v  v  ^ 


T=N.{(trK  2)  K  -  K  fl  K>.N, 


(3.28) 


0  =  det  K. 


-  13  - 


4.  Pure  homogeneous  deformation 

(a)  Propagation  in  principal  direction 

We  now  suppose  that  the  deformation  to  which  the 
body  is  subjected  is  the  pure  homogeneous  deformation ,  the 
principal  directions  for  which  are  along  the  axes  of  the  ref¬ 
erence  system  x.  Then 

ci ,  =  o  (i+j)  (4.1) 


and  it  follows  from  (2.3)  that 


k  =  0,  oKj  =  0  (i=}=j)  •  (4.2) 

The  principal  waves  are  waves  for  which  the  directions 
of  propagation  are  along  the  principal  directions  of  strain, 
i.e.,  the  waves  for  which 


n. 

n. 


or  6 .  , . 
i3 


(4.3) 


We  consider  first  the  waves  propagated  along  the  x^-axis. 
Then,  introducing  (4.2)  into  (3.16),  we  obtain 


kll  "  S%2' 


,k22  -  s2wi;i 


=  0, 


(4.4) 


whence 


s2  =  kU/“22  °r  s2  =  k22/Mll- 


(4.5) 


14 


We  assume  that  these  quantities  are  positive  and  consider 
the  waves  corresponding  to  the  positive  square  roots,  i.e., 
the  waves  travelling  in  the  positive  direction  of  the  x^-axis. 
We  employ  the  notation 


S13  ^kll//W22)  ^  '  S2 3  ~  (k22//°Jll) 


(4.6) 


We  note  from  (3.15)  that,  for  the  wave  for  which  s  =  s 

d_  =  b.,  =0  and,  for  the  wave  for  which  s  =  s00,  d.  =  b_  =  0. 

d  -L  d  5  JL  d 

Thus,  the  former  wave  is  polarized  with  d  and  b  in  the  x^ 
and  x_  directions  respectively  and  the  latter  with  d  and  b 
in  the  xg  and  x.^  directions  respectively.  It  follows  from 
(4.2)  that  for  these  waves  e  is  polarized  in  the  same  direction 
as  d  and  h  in  the  same  direction  as  b. 

More  generally,  we  adopt  the  notation  that  s.  .  (i=fj) 
is  the  slowness  for  the  principal  wave  whose  direction  of 
propagation  is  along  the  x^-axis  and  which  is  polarized  with 
its  electric  displacement  field  in  the  x. -direction.  Then, 
analogously  with  (4.6) ,  we  have  the  further  relations 


s32  (k33A>n)  /  s12  ' 

S21  ~  ^/i22^3'3)  t  '  s3i  =  'k33/’w22)  ^  ' 


(4.7) 


It  follows  from  (4.6)  and  (4.7)  that 


S13S32S21  S23S12S31 


(4.8) 


and  this  relation  is  valid  for  any  constitutive  equations  of 
the  form  (2.8)  with  k  and  w  given  by  (2.6),  k  ,  k  ,  k_  and 

.  ~  ®  ip  ^  p 

wo,  w  ,  w2  being  arbitrary  functions  of  tr  c,  tr  c  ,  tr  c  . 
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In  the  case  when  to.  =  to.  =  o  and  to  is  constant, 

1  2  o 

i.e.  the  specific  reluctance  is  independent  of  deformation, 

we  have  to....  =  to-n  =  to_.  =  to  and  it  follows  from  (4.6)  and 
11  22  33  0 

(4.7)  that 


S13  S12 '  S21  S23'  S32  S31* 


(4.9) 


If  on  the  other  hand  k,  =  k_  =  0  and  k  is  constant,  i.e. 
the  dielectric  constant  is  independent  of  deformation,  we  have 
ki;L  =  kg2  =  k^3  =  kQ  and  it  follows  from  (4.6)  and  (4.7) 
that 

S23  =  S32 '  S13  =  S31'  S21  =  S12* 


(4.10) 


We  shall  now  consider  the  somewhat  more  general 
case  when  the  direction  of  propagation  n  is  in  the  plane 
formed  by  two  of  the  principal  directions  of  strain.  Choosing 
the  coordinate  system  x  with  the  x2~axis  perpendicular  to  this 
plane,  we  have 


ni  =  (ni'°'n3> ' 


(4.11) 


and  the.  Finger  strain  components  are,  as  before,  c^,  with 

c.  =  0  (i=j=j) .  We  choose  a  new  coordinate  system  x  with  the 
^  j 

axis  x^  parallel  to  n  and  &2  coinciding  with  x^.  Then 


x.  =  a. ,x, , 
i  ij  y 


(4.12) 


where  a^  is  given  by 


ij 


n3,  0,  -nj 


0,  1,  0 


nx,  0,  n 


3  ii 


(4.13) 


The  components  c^  of  the  Finger  strain  tensor  in  the  system 
x  are  given  by 


i  Clln3  +  C33nl' 


0  ,  (c11-c33)n1n3 


^C11~C33^  nln3' 


'22  ' 


0  '  Cllni  +  C33n3 


(4.15) 


I  p  2  o  2 
!  Clln3  +  C33nl' 


0  '  (c^ii"c33)nin: 


CikCk j 


|  (c^1-c|,)n1n3. 


'22  ' 


2  2  2  2 
'  Cllnl  +  C33n3 


Referred  to  the  system  x,  the  dielectric  constant  and  inverse 
magnetic  permeability  matrices,  k  and  w  respectively,  are,  by 

•V 

analogy  with  the  relations  (2. 3), given  by 


k.  .  =  k  6.  .  +  k,c.  .  +  k_c.,  c,  . 
ij  o  xj  1  ij  2  lk  k j 


(4.16) 


u .  .  =  w  6  .  ,  U.  C .  .  +  U)„  c . ,  c.  .  , 
ij  o  ij  1  ij  2  ik  k j 


where  k  ,  k..  ,  . . . ,  to  have  precisely  the  same  meanings  as 
in  (2.3).  Since  tr  c,  tr  c  ,  tr  c3  are  invariant  under 
orthogonal  transformations,  k  ,  . ..,  w  may  be  expressed 

0  o  <- 

as  functions  of  tr  c,  tr  c“,  tr  c  of  the  same  forms  as  they 

~  p  ~  ^ 

are  functions  of  tr  c,  tr  c  ,  tr  c  .  Introducing  (4.15)  into 
(4.16) ,  we  see  that 

£*j  -  ■  5ij  ■  <5-1»  t j  ■  »• 

(ij  =  12,  21,  23,  32) .  (4.17) 

Introducing  (4.17)  into  the  secular  equation  (3.16) 
with  k  and  to  replaced  by  k  and  u>  respectively,  we  obtain 


s2  =  s23  (say)  =  ki;L/(c22 


(4.18) 


s2  =  s23  (say)  =  k22/w11. 


Again  assuming  k11/w22  and  k22/Si;L  are  both  positive,  we  see 
from  (3.15)  ,  with  k,  w,  d  and  b  replaced  by  k,  £5,  d  and  b 
respectively,  that  for  the  wave  for  which  s  =  s13,  dg  =  b1  = 
and  for  the  wave  for  which  s  =  sL„,  d.  =  b„  =  0.. 


5.  Simple  shear 


We  now  consider  the  case  when  the  deformation  is  a 
simple  shear  of  amount  K,  described  in  the  rectangular  cartesian 
coordinate  system  x  by 


=  X.^,  Xg  =  Xg,  x^  =  +  kX^  -  (5.1) 

We  consider  a  wave  propagated  in  the  positive  direction  of  the 
x^-axis,  so  that 


n.  -  6 . 
i  i3 


(5.2) 


From  (5.1)  and  (2.4),  we  obtain 


V 

o 

o 

K2/  0,  «3  j 

: 

o 

o 

o 

'  CikCk j 

o 

o 

o 

K,  0,  K2 

1 

i 

|  K3,  0  ,  K2  +  Ki| 
i  | 

Introducing  (5.3)  into  (2.3)  we  obtain 


(5.3) 


I  k  +k0K  , 
!  o  2 


kij  = 

1 


(kx+k2K  )  , 


,  K(k1+k2<  ) 


,  ko+K2(k1+k2+k2ic2)  | 


(5.4) 


W  +W0K  ,  0 

o  2 


,  K  (w^+WglC2)  I 


k(o>  +ft)  K  )  , 


2  2 
,  U)  +K  (U)-+C0o+W„K  ) 
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We  recall  that  k  ,  k.  ,  k_ ,  co  ,  ui, ,  w.  are  functions  of 
2  o'  1'  2'  o'  1  2 

tr  c,  tr  c  ,  tr  c5  winch  are  now  given  by 

tr  c  =  k2,  tr  c2  =  2k2  +  k*1,  tr  7  =  3k1*  +  k^.  (5.5) 


Introducing  (5.4)  into  (3.16),  we  obtain 


=  jiW2)/“o  1 


1/2  f~ 271  1/2 

or  _ kQ/ (u)o+032k  ^)J 


(5.6) 


From  (3.15)  we  see  that  the  first  of  these  values  of  s  leads 


d2  =  bl  =  0 


(5.7) 


and  the  second  leads  to 


dl  b2  =  °* 


(5.8) 


Thus,  in  the  case  provided  by  (5.6) 1,  d  and  b  are  linearly 
polarized  in  the  x-j^  and  x2  directions  respectively  and  in  the 
case  provided  by  (5.6) 2,  they  are  linearly  polarized  in  the 
x2  and  x1  directions  respectively.  The  corresponding  express¬ 
ions  for  e  and  h  can  be  obtained  by  introducing  (5.4)  and 
(3.14)  into  (2.8).  We  obtain,  corresponding  to  (5.7), 


e  =  const. 


k  +k2 (k  ■ 
o  1 


+k2+k2<2) ,  0,-k  (k1+k2<2)  , 


(5.9) 


h  =  const. 


[o,  X,  0_ 


21  - 


and,  corresponding  to  (5.8),  we  obtain 


e 


const. 


•J 


and 


(5.10) 


h  =  const. 


r— 

L* 


2 

o  2 


0, 


K  (0>1+t02K2) 


] 


22 


6.  Shear  in  two  directions 


We  now  consider  the  case  in  which  the  direction  of 
propagation  of  the  electromagnetic  wave  is  along  the  x^-axis, 
but  the  direction  of  shear  is  in  the  x2x^  plane,  the  deform¬ 
ation  being  described  by 

X1  ~  Xi>  x2  =  X2  +  ^xi'  X3  =  X3  +  kXi*  (6.1) 
Introducing  (6.1)  into  (2.4),  we  obtain 


o. 

X  , 

K 

- 

x. 

X2, 

Xk 

K, 

Xk, 

2 

K 

and  (6.2) 


C  .  ,  C,  . 
ik  kj 

= 

2  2 

X  +K  ,  0  ,  0 

0  ,  X2  ,  Xk 

i 

+  (X2+k2) 

0  ,  X  ,  K  jj 

i 

1 

X  ,  X2  ,  Xk i 

1 

OJ 

-< 

o 

2 

k  r  Xk  ,  k 

It  follows  from  (6.2)  and  (2.3)  that 

kQ  +  k2(X2+x2),  k  X  +  k2A(X2+K2)  i 

i 

k1X  +  kgX (X2+k2) ,  kQ  +  kxX2  +  k2X2(l+X2+K2) I 


and 


(6.3) 
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w 


aB 


= 


Wo  +  (02^2+k2^'  +  “p* (A2+k2) 


oj^X  +  u^X  ( X2+k2)  ,  wQ  +  a^X2  +  w2X2  (1+X2+k2) 


2  3 

k  ,  kn ,  . ..,  u  are  functions  of  tr  c,  tr  c  and  tr  c  ,  which 

U  Ji  £  <V  -w  <v 

are  given  by 


tr  c  =  X2+k2  ,  tr  c2  =  2(X2+k2)  +  (X2+k2)2  , 


tr  c3  =  3  (X2+ic2) 2  +  (X2+tc2)3. 


(6.4) 


Since  the  wave  considered  is  propagated  along  the 
x3  axis,  it  follows  that  the  relations  (3.14)  and  (3.15)  are 
satisfied.  We  have,  therefore. 


d3  =  b3  =  0, 


and 

(5ars2Aaeld8  “  0 

f 

(6.5) 

where 

A11  = 

“22  1 11  “21 <s  *  21' 

A12=  U)22{~  *  12 

W21(~  ) 22 
(6.6) 

A21  * 

“ll(~  ^  21  “  W12{~  *11' 

A22=  Wll{~  *  22 

W12  *  12 

From  (6.5) 2,  we  see  that 

s  is  given  by 

{ae  -  5  Aa(S 


i.e, 


=  0, 


(6.7) 


2  (A11+A22)  *  2  ^All"A22)  +  **A12A2l[  / 

(6.8) 
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=  T0  +  t2(X2+K2), 

22  =  To  +  Tlx2  +  t2X2(1+X2+k2)/  (6.14) 

1}12  =  (^-1)21  =  xiX  +t2X(X2+<2)  . 

Introducing  (6.14)  and  (6.3)  into  (6.6),  we 

obtain 


{WoVV2k2)  *  A2(wiTo+a,2To“WlTl+WoT2) 

+  A2K2(«2T0+«2T2-«2T1)  +  ^(T^Tg-Ti), 

(woTo+(,)2Tok2)  +  x2(“oTl+“oT2~(1,lTl+0''2To) 

+  X2k^wot2+w2t2-Wit2)  +  ^1|(woT2+a)2T2-0)1T2)  , 


=  X(a)oTra)iTo)  + 


(6.15) 


+  X  ( Vl^lV^oV^V  ' 


A21  "  X(“oTraJiTo)  +  X,c2(woT2"W2To+U,2Tl'a,lT2) 

+  ^3(w0Va,2Va,2VWlT2)‘ 

For  brevity  we  introduce  the  notation 


A11  = 

an  +  bnx2  + 

A12  = 

ai2X  +  b12X3, 

A21 

a21X  +  b21X  3 ' 

> 

ro 

ro 

II 

a22  +  b22X  +  C22X^ ' 

(6.16) 


where,  comparing  (6.16)  and  (6.15), 
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a 

b 


11 

11 


=  Vo  +  woV  ' 


=  (»1T0+»2T0-«1T1+«0T2)  +  K  (»2T0+«2T2-«2T1) , 


etc. 


(6.17) 


Introducing  (6,16)  into  (6.8)  ,  denoting  by  s^  and 
s2  the  values  of  s  obtained  by  taking  the  positive  and  negative 
square  roots  respectively  in  (6.8),  and  expanding  the  expres¬ 
sions  for  s^  and  s^  as  power  series  inX,  we  obtain 


expression 


=  a  +  b  X  +  a 


a22  +  ^22^  "  a 


Introducing  (6.16) 
obtained  as  a  power 


12a21 {ail"a22)  X  + 


- 1  2 

12a21(ail"a22)  X  + 


(6.18) 


into  (6.9)  and  expanding  the 
series  in  X ,  we  obtain 


(<Vdi} 

Xa21(ail"a22) 

s=sl 

(Vd2} 

Xai2 (a22~ail) 

s=s2 

(6.19) 


27 


7 .  Application  to  time-dependent  materials 

We  now  consider  that  the  dielectric  constant  matrix 
k  and  specific  reluctance  matrix  w  depend  not  only  on  the 
deformation  gradients  existing  at  the  instant  of  measurement, 
but  on  the  whole  history  of  the  deformation  gradients  in  the 
particle  considered  up  to  this  time.  This  means  that  k  and  w 
are  matrix  functionals  of  the  history  of  the  deformation  grad¬ 
ients.  However,  if  we  restrict  the  deformations  considered  to 
ones  in  which  the  body  is  taken  from  the  undeformed  state  to  a 
certain  state  of  deformation  at  some  instant  of  time  and  then 
held  there,  and  we  further  make  appropriate  assumptions  regard¬ 
ing  the  nature  of  the  functional  dependence  of  k  and  oo  on 
the  deformation  gradient  history  and  the  path  by  which  the 
material  is  taken  from  its  undeformed  state  to  the  steady  state 
of  deformations,  we  can  still  write  the  constitutive  equations 
for  d  and  h  in  the  forms  (2.5)  where  k  and  w  are  now  functions 
of  the  steady  state  deformation  gradients  and  the  time  which 
has  elapsed  since  these  were  produced  in  the  material.  All 
the  results  obtained  in  the  paper  then  follow  with  the  proviso 

that  (aQ ,  to1 ,  (o5  and  kQ,  k^  ,  k2  depend  on  this  time  as  well  as 

2  ^ 

on  the  strain  invariants  tr  c,  tr  c  and  tr  c  . 

This  parallels  the  application  of  results  in  finite 
elasticity  theory  to  problems  involving  stress  relaxation  in 
viscoelastic  solids  held  at  constant  deformation  [3,4). 
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8.  Appendix 

In  this  section,  we  outline  the  computation  which 
yields  the  secular  equation  (3.12).  With  (3.11),  we  have 


det  (k..+s  e.  e.  n  n  to  )  =  0, 

ij  ipq  jrs  p  s  qr 


(8.1) 


Let 


a .  =  sc.  e .  n  n  to 

ij  ipq.  jrs  p  s  qr 


(8.2) 


With  (8.2),  we  may  write  (8.1)  as 


6  det  (k. ,+a.  .)  =  e.  e,  (k.  .+a.  .) (k  +ct  )  (k  +a  ) 
lj  ij'  ipq  jrs'  lj  ij/v  pr  pr'  '  qs  qs' 

(8.3) 

=  (k.  .k  k  +3k.  ,k  a  +  3k..a  a  +a.  ,a  a  )  =  0, 

ipq  jrs  ij  pr  qs  ij  pr  qs  lj  pr  qs  lj  pr  qs 

We  then  have 


e.  e  k. .k  k  =6  det  k. 
ipq  jrs  ij  pr  qs 


(8.4) 


With  (8.2),  we  obtain 


3e.  e  . 
ipq  jrs 


kijkpraqs  =  2  [j}-(kuk).n  "  (n.k.n.)(tr  kto)j  , 


3£ipqejrskijW  ’  3s\n.k.n)  [  2n.<j2  .n+ (tr 


to) 


-  tr  to^"  -  2  (tr  to)  (n 


.  to.n)J  , 


(8.5) 


(8.6) 


and 


„a.  .a  a  =6  det  a..  = 
ipq  jrs  ij  pr  qs  ij 


0. 


(8.7) 
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With  (6.11) ,  we  note  that  the  expression  appearing  within  the 
brackets  in  (8.6)  may  be  written  as 

[  2  2  2  “!  -1 
2w  -(tr  u)d)+(tr  to)  I-(tr  w  )I  .n  =  2det  w(n.w  .n)  . 

(8.8) 


Substituting  (8.4),  . ..,  (8.8)  into  (8.3)  we  obtain 

det  (k.  +a.  .)  =  ^  -  ipS2  +  0  =  0  (8.9) 

i  j  i  j 

where 

<p  =  (n.k.n)  (n.ai_1.n)  det  a), 

ip  =  (n.k.n)tr  kw  -  n.kwk.n,  (8.10) 

0  =  det  k. 


Equations  (8.9)  and  (8.10)  yield  the  result  (3.12)  and  (3.13). 
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